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Abstract 

This paper presents techniques for digital image watermarking based on 
eigenvalue decomposition of normal matrices. The introduced methods are 
convenient and self-explanatory, achieve satisfactory results, as well as require 
less and easy computations compared to some current methods. Through the 
proposed methods, host images and watermarks are transformed to the space 
of normal matrices, and the properties of spectral decompositions are dealt 
with to obtain watermarked images. Watermark extraction is carried out 
via a procedure similar to embedding. Experimental results are provided to 
illustrate the reliability and robustness of the methods. 

Keywords: Digital image watermarking, Extracting, Embedding, Normal ma¬ 
trix, Eigenvalue. 


1 Introduction 

Nowadays, everyone has an easy access to digital images via the Internet and other 
means of communication. In addition, powerful software and hardware tools are 
available to edit any kind of digital multimedia. So, digital images may become a 
problematic issue if it have significant authentications and information. To cope 
with this problem, digital image watermarking can be employed. Digital image wa¬ 
termarking is a procedure of embedding additional information called watermark 
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into an image, called host image, which preserves the perceptual qualities of the 
original host image. When used normally, a watermark, carrying important infor¬ 
mation, is hidden and can be visible only by pursuing some particular processes. 
The watermark should be detectable or extractable from the host image by owner 
identihcation or integrity verihcation of the watermarked image [12]. 

Watermarking has a long list of applications in such areas as protecting copyright 
ownership of digital media, verifying the authenticity of a bill or the trademark of a 
paper manufacturer, making the information imperceptible, keeping the information 
secret, and so on. Digital image watermarking has received a great deal of attention 
as much as to be practiced in many other important applications. References for 
theory and practice of this method are HQllIIlill] and [13], to name but a few. 

For the subject to be fully grasped, presenting certain basic dehnitions seems 
necessary. A gray scale image can be demonstrated as a matrix the magnitude of 
each of whose components is called gray level or intensity of the image at that point. 
A given image is called a digital image if all of its elements are hnite and discrete 
quantities; in this case, these elements are called pixels. Digital Image Processing 
(DIP) concerns using digital computers to preform processing on digital images. 
DIP receives a digital image as an input, processes efficient algorithms on it by a 
computer and gives an image as an output. So, digital image watermarking can be 
construed as a DIP. 

The theory of linear algebra provides helpful instruments for digital image wa¬ 
termarking as well as a variety of other helds of digital image processing. One of 
the most useful tools for digital image watermarking is employing Singular Value 
Decomposition (SVD) [3lll]. 

the matrix A G can be written in the form of USV*, where U G 

and V G are unitary matrices, i.e., U*U = Im,V*V = In, where * denotes 

complex conjugate transpose. The matrix S' is an m x n diagonal matrix such 
that its nonnegative entries are ordered in a non-increasing order (see for example. 
Theorem 7.3.5 of [S]). Motivated by the influences and traits of singular values of A 
in digital image watermarking, and regarding the point that the singular values of A 
are the positive square roots of the eigenvalues of matrices AA* and A* A, this paper 
concerns itself with the eigenvalue of the normal matrices A + A* and A — A* in 
order to obtain techniques for digital image watermarking which are more effective, 
yield better results and need fewer computations. 

In the next section, some dehnitions and concepts are briehy presented about 
normal matrices. Section [3] consists of two subsections in which the proposed image 
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watermarking methods are explained. In Section 01 the validity rates of the pre¬ 
sented watermarking schemes are investigated, and their efficiencies are compared 
by experimental results. In addition, the robustness of watermarking is evaluated 
against some common attacks and criticisms. 


2 Normal matrices 


In this section, definitions and some properties of normal matrices are reviewed. 
See [21 [8] and references therein as the suggested sources on a list of conditions 
on normal matrices. The proposed method based on these presented properties is 
explained in the next section. 

A matrix M G C"^"' is called normal if M*M = MM*. Assume that M is an 
n-square normal matrix, then there exists an orthonormal basis of consisting 
of eigenvectors of M, and M is unitarily diagonalizable. Let the scalars Ai,..., A„, 
counted according to multiplicity, be eigenvalues of the normal matrix M and let 
Ml,... be its corresponding orthonormal eigenvectors. Then, the matrix M can 
be factored as 


M = UAU*, A = diag(Ai,..., A„), U = [ui,... ,Un\, 

where the matrix U satishes UU* = In- Without loss of generality, assume that 
eigenvalues are ordered in a non-ascending order of magnitude, i.e., |Ai| > IA 2 I ... > 

|A„|. 


Note that, if all the elements of the matrix M are real, then M* = M^, where 
denotes the transpose of the matrix M. A square matrix M is called symmetric 
if M = M'^ and called skew-symmetric if M = —M^. It is easy to see that sym¬ 
metric and skew-symmetric matrices are normal. Also, all the eigenvalues of a real 
symmetric matrix are real, and all the eigenvalues of a real skew-symmetric matrix 
are purely imaginary. A general square matrix M satishes M = B C, where the 

M + M^ 

symmetric matrix B = - is called symmetric part of M and, analogously. 


the skew-symmetric matrix C = 


M - M^ 


is called skew-symmetric part of M. Con¬ 


sequently, every square matrix can be written as the sum of two normal matrices: 
a symmetric matrix and a skew-symmetric one. This point is specially used in the 
proposed watermarking schemes. 
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3 Digital image watermarking method 


This section includes two subsections in which we present methods for digital image 
watermarking using normal matrices. To do this, we transform both the matrix 
representing the host image and the one representing the watermark into the space 
of normal matrices. Next, we benefit from the properties of its eigenvalue decom¬ 
positions and embed the watermark into the host image. Finally, returning to our 
original space, we can construct the watermarked image. Also, we may extract the 
watermark by following almost the same processes. 

Let X and W be two square matrices representing the host image and the water¬ 
mark image respectively. Two distinct methods are considered. First, we deal with 
the symmetric parts of X and W to obtain a digital image watermarking scheme. 
Next, we explain another technique using both symmetric and skew-symmetric parts 
of the matrices X and W. It is noticeable that hnding the eigenvalues and eigen¬ 
vectors of a matrix needs less computation than finding its singular values and 
singular vectors. Furthermore, the eigenvalues and eigenvectors of a normal (espe¬ 
cially symmetric or skew-symmetric) matrix can be calculated by explicit formulas 
and, therefore, may yet again require less computation m- 

3.1 Digital image watermarking method using symmetric 
parts of X and W 

In this subsection, we describe a digital image watermarking scheme based on the 
eigenvalue decomposition of the symmetric parts of the matrices X and W. This 
procedure can be carried out similarly for the case of the skew-symmetric parts of 
the two matrices X and W. Notice that, for the same host image and watermark, 
results obtained by these two techniques (using symmetric or skew-symmetric parts) 
may be different. 

Assume that X is an n x n matrix and IF is an m x m. Let Bx and Bw be the 
symmetric parts of the two matrices X and IF, respectively. Therefore, the normal 
matrices Bx and Bw can be factored as follows 

Bx ^ ^Bx * * * ? 

Bw = Ubw^Bw^Bw-- = diag(AB^4,..., \Bw,n)- 

Let a be a given scaling factor. The parameter a regulates the strength of the 
embedding watermark into the host image. Now, embed the symmetric part of the 
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watermark into the symmetric part of the host image by modifying the eigenvalues 
of Bx and then obtain the symmetric part of the watermarked image Y as detailed 
below: 

By = Ubx-^yU'^^, Ay = Abx + OiKb^- (1) 

Hereinafter, reconstruction of the watermarked image Y from its symmetric part 
is rendered. To do this, all the elements that are located above the main diagonal of 
the matrix X are required. Since Y can be construed as an acceptable approximation 
of X, let us consider all the elements above of the main diagonal of Y as the elements 
of X. Obviously, the elements located below the main diagonal of the matrix Y^ 
are also determined. In addition, it is concluded from the fact 2By = Y + that 
the elements which are stated below the main diagonal of Y can be obtained by 
subtracting the elements stated below the main diagonal of 2By from the elements 
of Y'^ . See the above equation where ” /” and ” x ” denote the given and unknown 
entries, respectively. 


2 By — 

X 

/ ■ 

+ 

X 

X 


X 

X 


/ 

X 


Y yr 


Finally, it is clear that the elements on the main diagonal of Y are the main 
diagonal entries of By- It should be noted that by this procedure, the watermark 
is embedded only into the elements that are placed below the main diagonal of X, 
and its other elements are kept unchanged. 

Furthermore, to reconstruct the watermarked image Y, we may reserve the 
elements located below the main diagonal of X instead of the elements above the 
main diagonal and then follow a procedure similar to what is performed newly. 
Indeed, we can partition X into several segments and reserve it provided that Y can 
be reconstructed. This may be useful specially when unchanging some partitions of 
the host image is desirable during the digital image watermarking. 

Now, as the next necessary step, we must be able to extract the watermark from 
the watermarked image. Our scheme for watermark extraction implicates obtaining 
the symmetric matrix Bw and then reconstructing the watermark image W from 
it. In order to compute the matrix Bw, the diagonal matrix Ab^ and the unitary 
matrix Ub^ are required. Keeping ([T]) in mind, the following relations describe our 
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proposed extraction technique for finding By/ 


= 


Ay — Kbx 
a ’ 


Bw - Ubw^BwU*b^- 


Anew, reconstructing the W needs the elements which are stated above (or be¬ 
low) the main diagonal of W and can be performed similar to what was mentioned 
in the case of reconstructing the watermarked image Y. Let us once again remark 
that, if some segments of the watermark have more signihcance, then we can par¬ 
tition the watermark such that the important segments be preserved provided that 
reconstruction of the watermark is done. 


3.2 Digital image watermarking method using both sym¬ 
metric and skew-symmetric parts of X and W 

In the previous subsection, a digital image watermarking scheme which embedded a 
watermark image into almost half of a host image was introduced. That technique 
kept more than half of the elements of the host image unchanged, which might lead 
to some usages and advantages and, in particular, was a suitable method for small 
enough values of the scaling factor a. However, embedding all of the watermark into 
half of the original image may cause the watermarked image to lose its reliability as 
compared to some other watermarking schemes for large enough scaling factors. To 
cope with this defect, we introduce a new method concerning both symmetric and 
skew-symmetric parts of the watermark and the host image in order to embed the 
watermark in the host image completely. This new watermarking technique provides 
for a remarkably high reliability. 

As already mentioned, every matrix is equal to the sum of its symmetric and 
skew-symmetric parts. In the proposed technique, therefore, a given scaling factor is 
shared equally between the symmetric and skew-symmetric parts of the host image. 
Keeping in mind the definitions of Bx and Byy, Cx and Cw may be employed to 
denote the skew-symmetric parts of X and W, respectively. The matrices Cx and 
Cy/ may be written as follows: 

Cx = Ucx^Cx^Cx-' ^Cx = diag(Ac^4,..., \cx,n), 

^Cw^Cw^Cw ^ . . . , ^Cw ' 
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Inspired by the method described in the previous subsection, for scaling factor 
a, set 

Cx = a = Acx + 

~ ~ ~ Of 

Bx = UbxAUb^, a = Abx + 2^Bw- 

Therefore, the watermarked image X will be A” = Bx + Cx- 

Henceforth, extracting the watermark W will be considered to be from the wa¬ 
termarked image X. The extraction process needs two diagonal matrices, i.e., Abx 
and Acx, and two unitary matrices,i.e., Ub^^t and the presented water¬ 

mark extraction scheme, at first, the matrices A^^ and Acy^r are computed as the 
estimations of Ab^ and Ac^, respectively. Next, the extracted watermark W will 
be obtained by approximating the two matrices By/ and Cw Assume now, X is 
an acceptable approximation of the original image X. Form the symmetric and 
skew-symmetric parts of X, namely Bx and Cx, in the order already mentioned, 
the eigenvalue decompositions can be computed as shown below: 

X + X'^ 

Bx = - 2 -’ Bx = Ubx-^BxUbx^ 

= Cx = Uc,kc,Uh,. 

Now, the diagonal matrices A^^ and Aq^^ are obtained from the following rela¬ 
tions: 

X _ ^Bx T _ ^Cx ^Cx 

a a 

Also if we consider Bw and Cy/ estimations of By/ and Cw, respectively, then 
Bw = Ubv^t^BwBbw' Cw = Uc^Acy^Uc^. 

Finally, the extracted watermark W will be W = Bw + Cw- 

We should note that none of the four required matrices in the extraction process 
provides information about the host image and the watermark image. This point 
improves the safety of the presented watermarking method. 
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4 Experimental results 


In this section, the validity and the influence of the watermarking methods de¬ 
scribed in Section [3] are examined by experiments. Moreover, let us note that the 
ideas presented in this paper, can be used immediately to establish a block-based 
watermarking technique. A block-based watermarking scheme concerns dividing 
both the watermark and the host image into non-overlapping blocks and applying a 
watermarking method to each block. The block-based SVD watermarking method 
is especially well known (e.g., see Section 4 of [T]). 

In our experiments, the watermarking method explained in Subsection 13.11 is 
considered for using either the symmetric or skew-symmetric parts of the host image 
and the watermark. Also, the watermarking scheme demonstrated in Subsection 13.21 
has been tested to prove its efficiency. 

The Peak Signal to Noise Ratio (PSNR) is calculated to measure the quality of 
the watermarked image. PSNR for gray scale images of size M x N, whose pixels 
are represented with 8 bits, is computed by the following relations: 

255 ^ 1 

PSNR = 10 logjf®; MSE ^- Ntjf, 

ij 

where Xij and Ajj denote the elements of the original and the watermarked images, 
respectively. In the above relation, MSE denotes the Mean Square Error between 
the host image and the watermarked image pixels. In the performed experiments, 
we considered four host 512 x 512 gray scale images (a) Lena, (b) Cameraman, 
(c) Baboon and (d) Peppers presented in Fig [T] Also, our watermark image was 
Yazd University emblem of size 128 x 128, as displayed in Fig [2l The results are 
shown in Tables [H El |3] and IH for some values of the scaling factor a, for the host 
images Lena, Cameraman, Baboon and Peppers and the watermark image Yazd 
university emblem, respectively. In these tables. Method 7^1 and Method 7^2 stand 
for using the symmetric and the skew-symmetric parts of the host and watermark 
images, in the image watermarking schemes, respectively. In addition. Method 7^3 
denominates the watermarking technique described in Subsection 13.21 The results 
obtained by these techniques are compared to the results achieved by the digital 
image watermarking method using SVD which is denoted by Method 7^4 in our 
tables. 

Furthermore, the robustness of our watermarking schemes was evaluated in the 
presence of different types of attacks. The Bit Error Rate (BER) was used to mea¬ 
sure the robustness between the original watermark and the extracted watermark. 









BER is the percentage of bits with errors divided by the total number of bits that 
have been processed. We should note that the proposed methods described in Sub¬ 
section 13.11 are highly robust, and the BER results for those two methods almost 
zero under many types of attacks. Hence, the results corresponding to the other 
watermarking method, which was denoted by Method ^3, are only presented and 
compared with results obtained by Method 7^4 to illustrate the robustness of our 
technique. The attacks used for examining the robustness of the method are median 
filtering, Gaussian low pass hltering, averaging filter, wiener hltering, JPEG com¬ 
pression, additive Gaussian white noise (AGWN), salt and pepper noise, resizing 
and adjust image intensity values. Table [5] shows the BER results for the host im¬ 
age Lena and the watermark image Yazd University emblem under various attacks, 
where the window size is 3 x 3. Moreover, Fig 12] shows the extracted watermark 
without any attack as well as after applying some sample attacks. 


Table 1; PSNR results for Lena and Yazd University emblem 


a 

Method#! 

Method#2 

Method#3 

Method #4 

0.2000 

72.1125 

84.1944 

74.8323 

74.8500 

0.4000 

66.0919 

78.1738 

68.8117 

68.8294 

0.6000 

62.5701 

74.6520 

65.2899 

65.3076 

0.8000 

60.0713 

72.1532 

62.7911 

62.8088 

1.0000 

58.1331 

70.2150 

60.8529 

60.8706 

1.2000 

56.5495 

68.6314 

59.2693 

59.2870 

1.4000 

55.2106 

67.2925 

57.9304 

57.9481 

1.6000 

54.0507 

66.1326 

56.7705 

56.7882 

1.8000 

53.0277 

65.1096 

55.7475 

55.7652 

2.0000 

52.1125 

64.1944 

54.8323 

54.8500 


5 Conclusion 

In this paper, the eigenvalue decomposition of symmetric and skew-symmetric ma¬ 
trices, as two important kinds of normal matrices, and their properties were used to 
obtain digital image watermarking schemes. The proposed methods were straight¬ 
forward and uncomplicated methods requiring clear and less computations compared 
to some exiting methods. This watermarking technique is the first one ever intro¬ 
duced and concerns the symmetric parts of a host image and a watermark, which 
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Fig 2: The Watermark Image 
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0 0 0 

(a) (b) (c) 

0 0 0 

'^! '-' 

(d) (e) (f) 

Fig 3: Extracted watermark from watermarked image (a) Without any attack (b) 
After median filtering attack (c) After Gaussian low pas filtering attack (d) After 
JPEG compression (Quality=50) attack (e) After salt and pepper noise (density of 
0.01) attack (f) After resizing (512-384-512) attack to the watermarked image. The 
used host image is Lena. 
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Table 2: PSNR results for Cameraman and Yazd University emblem 


a 

Method#! 

Method#2 

Method#3 

Method #4 

0.2000 

72.1120 

84.1944 

74.8323 

74.8500 

0.4000 

66.0914 

78.1738 

68.8117 

68.8294 

0.6000 

62.5696 

74.6520 

65.2899 

65.3076 

0.8000 

60.0708 

72.1532 

62.7911 

62.8088 

1.0000 

58.1326 

70.2150 

60.8529 

60.8706 

1.2000 

56.5490 

68.6314 

59.2693 

59.2870 

1.4000 

55.2101 

67.2925 

57.9304 

57.9481 

1.6000 

54.0502 

66.1326 

56.7705 

56.7882 

1.8000 

53.0272 

65.1096 

55.7475 

55.7652 

2.0000 

52.1120 

64.1944 

54.8323 

54.8500 


can be also considered for the case of skew-symmetric parts of aforesaid images. It 
has proved to be highly robust method against many types of attacks. This method 
could keep half (or some selected segments) of the image unchanging and this feature 
may be applied for some usages. However, this watermarking method could embed 
the watermark into just half of the original image, which caused the watermarked 
image to lose its reliability for large enough values of the scaling factor a. Hence, 
the second watermarking scheme using both symmetric and skew-symmetric parts 
of the original image and watermark was proposed. The experimental results illus¬ 
trate the high reliability and robustness of the methods. Finally, it is to be pointed 
out that these methods can be used to devise block-based watermarking techniques. 
Applying normal matrices to establish an image compression method will be the 
subject of our future research. 
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Table 5: BER results for the host image Lena and the watermark image Yazd 
University emblem under sample attacks 


Attack Method #3 Method #4 


Median hltering 

0.2319 

0.3479 

Gaussian low pas hltering 

0.0793 

0.2075 

Averaging hlter 

0.5249 

1.0498 

Wiener hltering 

0.2869 

0.3296 

JPEG compression (Quahty=70) 

0.1587 

0.0183 

JPEG compression (Quahty=50) 

0.2869 

0.0305 

AGWN 

0.3540 

0.8423 

Salt and pepper noise (density of 0.01) 

0.5005 

1.2939 

Resizing (512-384-512) 

0.0427 

0.1099 

Resizing (512-1024-512) 

0.0122 

0.0427 

Adjust image intensity values 

0.0427 

66.0889 
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